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We study an intersecting D-brane model which at low energies describes (l + l)-dimensional 
chiral fermions localized at defects on a stack of N c L>4-branes. Fermions at different 
defects interact via exchange of massless (4 + l)-dimensional fields. At weak coupling 
this interaction gives rise to the Gross-Neveu (GN) model and can be studied using field 
theoretic techniques. At strong coupling one can describe the system in terms of probe 
branes propagating in a curved background in string theory. The chiral symmetry is 
dynamically broken at zero temperature and is restored above a critical temperature T c 
which depends on the coupling. The phase transition at T c is first order at strong coupling 
and second order at weak coupling. 
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1. Introduction 

In this paper we continue the investigation of dynamical symmetry breaking in inter- 
secting D-brane models. The construction of these models starts with N c BPS Dp-branes, 
which we will refer to as color branes. The low energy theory on these branes is (p + 1)- 
dimensional U(N C ) super Yang-Mills (SYM) with sixteen supercharges. For p > 3, the 
Yang-Mills coupling g p +i scales like a positive power of length, and the SYM degrees of 
freedom become free in the infrared. 

To get non-trivial infrared dynamics we add additional D-branes, which we will refer 
to as flavor branes. The flavor branes intersect the color ones on subspaces of the (p + 1)- 
dimensional worldvolume of the latter. Open strings stretched between the color and flavor 
branes give rise to light fermions (and sometimes bosons as well) that are localized at the 
intersections. At low energies, fermions at different intersections interact via exchange 
of (p + l)-dimensional massless fields living on the color branes. In some cases, these 
interactions lead to non-trivial infrared effects. 

When the (p + l)-dimensional U(N C ) SYM theory is weakly coupled in the infrared, 
the effective interaction between fermions living at different defects becomes weaker as the 
distance between the defects, L, increases. The leading long distance dynamics can be 
studied in an approximation where one includes the effective four-Fermi interaction that 
arises from exchange of a single (p + 1) -dimensional massless field between two fermions, 
but neglects all the higher order exchange processes. Such processes are suppressed since 
the gauge coupling g p +i goes to zero at long distances. 

A nice thing about this approximation is that it leads to a theory that is tractable 
at large N c . Indeed, since the only dynamical degrees of freedom are the fermions, which 
transform in the vector representation of U(N C ), the low-energy dynamics is described by a 
vector model, and can typically be solved using familiar large N c field theoretic techniques. 

As L decreases, the interactions between the fermions and those between the fermions 
and the color degrees of freedom become stronger. For small L the dynamics can in some 
cases be analyzed using holography. In this limit the N c color branes can be replaced by 
their near-horizon geometry, and the interactions among the fermions are encoded in the 
dynamics of the flavor branes in the resulting curved spacetime. 

An example of an intersecting brane system of the kind described above is the DA — 
D8 — D8 system, which was studied from the present perspective in following some 
closely related earlier work |2|,[|. In that model the field theoretic analysis of the weak 
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coupling regime is complicated by the fact that the effective four-Fermi interaction due 
to single gluon exchange is long range - it is given by a non-local generalization of the 
Nambu-Jona-Lasinio model This leads to some interesting subtleties in the above 
considerations to which we hope to return elsewhere. 

The main purpose of this paper is to analyze an intersecting brane system which does 
not suffer from such subtleties. Our discussion is organized as follows. 

In section 2 we describe a brane configuration consisting of N c color D4-branes and 
Nf flavor D6 and _D6-branes in type IIA string theory, intersecting on a (1 + 1) -dimensional 
spacetime. The DA — D6 and DA — D6 intersections are separated by a distance L in the 
1R 3 along the fourbranes and transverse to the sixbranes. At the two intersections one 
finds chiral left and right-moving fermions, respectively. The low-energy dynamics of these 
fermions is governed by exchange of massless fields that live on the D4-branes, which are 
described by a five-dimensional gauge theory with 't Hooft coupling A (which has units of 
length) . The strength of the interaction among the fermions is determined by the value of 
the dimensionless parameter X/L. 

In section 3 we study the system at weak coupling, where one can (to leading order 
in X/L) restrict to the single gluon exchange approximation. The resulting four- Fermi 
interaction between the left and right-moving fermions is described by the Gross-Neveu 
model 0, with a specific UV cutoff. This model is known to break chiral symmetry, and 
to generate a mass scale. We rederive the chiral symmetry breaking from our perspective, 
and show that the dynamically generated mass of the fermions is much smaller than 1/L. 
This allows us to take a decoupling limit in which the intersecting brane model reduces 
precisely to the GN model. 

In section 4 we discuss the system at strong coupling, i.e. for A ^> L. In this regime 
the fermions and massless fields on the L>4-branes are strongly interacting, but there is still 
a weakly coupled description of the physics. It involves replacing the color D4-branes by 
their near-horizon geometry, and studying the flavor D6-branes as probes in this geometry. 
As in |IJ , there are two possible shapes that the flavor branes can take. One is the original 
separated parallel D6 and D6-branes. The other is a single connected brane which looks 
like the original D6 and L>6-branes connected by a wormhole whose width increases with 
X/L. We find that (as in [TJJ) the curved, connected brane has lower energy than the pair of 
separated straight ones. This implies that the chiral symmetry acting on the left and right- 
moving fermions is dynamically broken in this regime as well. The energy scale associated 
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with the breaking increases with the coupling X/L. We also comment on the spectrum of 
the model in the strong coupling regime and its relation to that of the Gross-Neveu model. 

In section 5 we discuss the finite temperature behavior of the model. In the weak 
coupling regime the chiral condensate decreases as the temperature increases and vanishes 
above a certain critical temperature T c . Thus, the chiral symmetry is restored in a second 
order phase transition at T = T c . At strong coupling, the (free) energy difference between 
the connected curved and separated straight branes decreases as the temperature increases, 
and beyond a certain critical value of the temperature the straight branes have lower free 
energy. The chiral symmetry restoration transition is strongly first order in this limit. 

We end in section 6 with conclusions and a discussion of some open problems. 

2. DA — D6 — D6 brane configuration and its low-energy dynamics 

The brane configuration that we will consider is depicted in figure 1. It contains the 
following D-branes: 

0123456789 
D4 : x x x x x (2.1) 
.D6, D6 : x x x x x x x 

N c color D4-branes stretched in (01234) are located at the origin in (56789). Nf flavor D6 
and L>6-branes are stretched in (0156789) and separated by a distance L in (234). Without 
loss of generality we can take the separation between the flavor branes and anti-branes to 
be in the x 4 direction. The color and flavor branes intersect on the (1 + l)-dimensional 
spacetime labeled by (01). We will be primarily interested in the low-energy dynamics at 
the intersection. 

The (9 + l)-dimensional Lorentz symmetry is broken by the branes ( |2.1|) as follows: 

SO(l, 9) -> SO(l, l)oi x SO(2) 23 x SO(5) 56 789 • (2.2) 

Each of the two intersections preserves a chiral supersymmetry. The D4 — D6 intersection 
preserves eight right-handed supercharges under SO(l, 1), while the D4 — D6 intersection 
preserves eight left-handed supercharges. The full system breaks all supersymmetry. 

We will study the system of figure 1 in the limit N c — > oo, g s — > 0, with g s N c and Nf 
held fixed and L>I S . In this limit the coupling of 6 — 6, 6 — 6 and 6 — 6 strings goes to 
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Fig. 1: The configuration of N c .D4- branes and Nf D6 and D6-branes which 
reduces at weak couphng to the Gross-Neveu modei. 

zero and they become non-dynamical external sources. The gauge symmetry on the D6 
and -D6-branes, 

U(N f ) L x U(N f ) R , (2.3) 

gives rise to a global symmetry of the (1 + 1) -dimensional theory at the intersection. 

The low-energy degrees of freedom are low lying 4 — 4 strings stretched between color 
branes, as well as 4 — 6 and 4 — 6 strings stretched between color branes and flavor branes 
and anti-branes, respectively. The low energy dynamics in the 4 — 4 sector is described 
by dimensional reduction of (9 + l)-dimensional N = 1 SYM with gauge group U[N C ) to 
4 + 1 dimensions. In the 4 — 6 sector one finds left-moving fermions qL transforming in the 
[Nf, 1) of the global symmetry (|2.3| ). 4 — 6 strings give right-moving fermions q R which 
transform as (1,7V/). Both qL and qn transform in the fundamental (N c ) representation 
of the U[N C ) gauge symmetry on the D4-branes. Note in particular that the U[Nf)i, 
symmetry in (|2.3|) acts only on the left- moving fermions, while U[Nf)n acts only on the 
right-moving ones. Thus, (|2.3| ) is a symmetry of the Lagrangian, but we will see later that 
it is broken to the diagonal U[Nf) by the dynamics. 

The 't Hooft coupling of the (4 + l)-dimensional gauge theory on the D4-branes, A, 
has units of length. The theory is weakly coupled for distances much larger than A and 
strongly coupled otherwise. In order to define it for all distance scales one must supply a 
UV completion. String theory does that by embedding it in a six-dimensional CFT, the 
(2, 0) theory. The (4 + 1) -dimensional gauge theory is the low-energy limit of the (2, 0) 
theory compactified on a circle whose radius is proportional to g s l s . 
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The fermions arising from 4 — 6 and 4 — 6 strings live on co dimension three defects in 
the (4 + l)-dimensional gauge theory. Their separation, L, introduces a second distance 
scale into the dynamics. One can think of the defects as probing the dynamics of the gauge 
theory (or, more properly, of the (2, 0) theory compactified on S 1 ) at the distance scale L. 

For L ^> A the theory on the fourbranes is weakly coupled. Naively, the (4 + 1)- 
dimensional fields that transform in the adjoint of SU(N C ) decouple in this limit, just like 
the corresponding modes on the sixbranes. This is not quite correct, since exchange of 
(4 + l)-dimensional massless modes leads to a weak attractive interaction between the left 
and right-handed fermions for any finite A. As we will see, this interaction leads to chiral 
symmetry breaking and has to be kept in the low energy description. 

The effective action that describes the dynamics in this limit contains the fermions 
Ql, Qr, and five-dimensional gauge fields on the D4-branes, Am, M = 0,1, 2, 3, 4, 

S = Jd 5 x -^Fl IN +d 3 (x-x 1 )q[a^(id^ + A^q L + 5 3 (x-x 2 )q R a fl (id fl + A fl )q R , 

95 (2-4) 
where the notation is as follows. The (4 + l)-dimensional worldvolume of the color D4- 
branes is parametrized by {x M } = x 1 , x). The flavor branes and anti-branes are placed 
at x = x\ and x 2 , respectively. As mentioned above, we can take X\, x 2 to lie on the x 4 
axis, with \x-\_ — x 2 \ = L. The indices \i in (|2.4| ) run over the range \i = 0, 1 or +, — . q^, q R 
are complex one component spinors. The corresponding one-dimensional Dirac matrices 
areo** = <y£, a M = St. 

We will define the 't Hooft coupling of the (4 + 1) -dimensional gauge theory on the 
color branes in terms of g$ (|2.4|) by0 

If- < 2 - s > 

As mentioned above, A is a length scale. In the next section we will use the Lagrangian 
( p.4j) to analyze the interaction between the fermions q^ and q R for A < L. In the opposite 
limit of large X/L we cannot use (pT4]), both because we have to add to it the other massless 
fields living on the D4-branes, and because in that regime the dynamics of these modes is 
no longer described by weakly coupled field theory in 4 + 1 dimensions. Nevertheless, one 
can still efficiently study the theory in this limit. This will be discussed in section 4. 



1 Here and below we often set a' = 1. 
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3. Weak coupling analysis 

In order to study the dynamics of the fermions qL, qR at weak coupling, it is convenient 
to integrate out the five-dimensional fields living on the color L>4-branes. To leading order 
in X/L this can be done by taking into account the interaction among the fermions due to 
single gluon exchange, using the action ( |2.4j ). The calculation is similar to that of |IJ] and 
leads to 



<S e s =i J d 2 x (qi^d^L + q R a^d^q R 
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4tt 2 



J d 2 xd 2 yG(x - y, L) (q ] L {x) ■ q R {y)^j (q R {y) • qh{x) 



(3.1) 



As in q\(x) ■ qniy) is a global color singlet (hence the dot) and transforms in the 
(Nf,Nf) of U(Nf)L x U(Nf) R . G(x,L) is the five-dimensional massless propagator over 
a distance x in the directions (01) and L in the directions (234), 

G(x,L)= 1 ■ (3.2) 
(x z + L z ) 2 

In addition, there are self interactions between fermions at a given intersection (the analog 
of eq. (3.2) in [p]), but these are small for distances much larger than A and can be 
neglected. We also Wick rotated to Euclidean space, which is convenient for analyzing the 
vacuum structure. 

Many properties of the solution of the interacting quantum field theory corresponding 



to the action (|3.1| ) can be anticipated as follows. The main complication in ( |3.1| ) is the non- 
local four Fermi interaction on the second line. If the function G(x — y, L) was proportional 
to 5 2 (x — y), ( |3.1| ) would reduce to the Gross-Neveu model 0, an asymptotically free local 
quantum field theory which is exactly solvable in the large N c limit. 

The actual Green function, ( |3.2|) can be thought of as a 5-function smeared over a 



region of size L. Using the fact that 



d 2 xG(x,L) = ^, (3.3) 
we can, at length scales large compared to L, replace ( |3.1|) by the local action 
Seff = J d 2 x iq[w^d^q L + iq^d^R + (ql( x ) " Qr( x )) (<1r( x ) ' Ql(x)^ (3.4) 
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with a UV cutoff A ~ -t. Comparing to we see that the Gross-Neveu couphng X gn is 
given in terms of our parameters by 

*gn — —j— ■ (3-5) 



In m it was found that at large N c the model fl3.4|) dynamically generates a fermion mass 
which goes like 

M F ~ fie A f " (3.6) 

where \i is an arbitrary renormalization scale and X gn the coupling at that scale. Since 
our model reduces in a particular limit (which will be made more precise below) to the 
Gross-Neveu model, it is natural to expect that it exhibits similar behavior. We next verify 
this by a more detailed analysis following closely that of [0. 

To solve the theory ( |3.1| ) for large N c it is convenient to introduce a bilocal field 
T(x, y) which is a singlet of global U(N C ) and transforms as (Nf, Nf) under the symmetry 
(p|). The action flOD can be rewritten as 



S eS = i J d 2 x (qicr^d^qL + q^a^d^R 

N c T(x,y)T(y,x) 



+ I d 2 xd 2 y 



+ T{y, x)q\(x) ■ q R {y) + T(x, y)q R {y) • ?l( 



x 



A G(x — y, L) 

(3.7) 

Integrating out T(x, y) using its equation of motion, 

T (x, y) = ^~ G ( X - y> L )qi( x ) ■ qr(v) > ( 3 - 8 ) 

one recovers the original action ([O]). To solve the large iV c theory one instead integrates 
out the fermions and obtains an effective action for T, which becomes classical in the large 
N c limit. As in Pj], using translation invariance of the vacuum, it is enough to compute 
this effective action for T(x, y) = T(\x — y\). Dividing by N c and the volume of spacetime 
we get the effective potential 

where T(k) is the Fourier transform of T(x). Varying ( |3.9| ) with respect to T(k) leads to 
the gap equation 

XJ G(x,L) k 2 + T(k)T(k) v ' 



The solution of (|3.10| ) gives the expectation value of the chiral condensate (q^ix) ■ qii(y)) 
via ( p.8|) . This equation has a trivial solution T = 0, but we will next see that there is a 
non-trivial solution as well. As shown in |JJ , any well behaved non-trivial solution of ( |3.10| ) 
has a negative value of V e g; hence the vacuum of the theory breaks chiral symmetry. 

To solve the gap equation for small X/L it is convenient to define the field f(x, y) by 

T(x : y) = ^-G(x-y,L)f(x,y) . (3.11) 

Using (|3.8| ) we see that the expectation value of f(x,y) is proportional to the condensate 

f{x -y) = {f{x, y)) = ^(ql(x) ■ q R {y)} . (3.12) 

In particular, the condensate at coincident points, {q L (x) ■ qn{x)) is proportional to /(0), 
which we will denote by m/. We will determine my as part of the solution of the gap 
equation. We will assume (and justify later) that f(x) is approximately constant on the 
scale L. Since the function G(x, L) fl3.2|) goes rapidly to zero for x > L (and can be thought 
of as a smeared 5-function) this implies that 

T{k) = — [ d 2 xG(x,L)f(x)e- lk - x ~ m /e " fcL . (3.13) 
2ir J 

Substituting (|3.11|) and ( |3.13| ) into the gap equation (|3.10| ) gives the Fourier transform of 

m = • (3-14) 

We can now determine the mass parameter m/ by evaluating /(0) in terms of its Fourier 
transform: 

d 2 k mfe~^ L 



m f = /(0) = \ gn J {2n)2k2 + m 2 e - 2m ■ (3-15) 
Dividing by my and estimating the momentum integral for m/L < 1 we find 

l~^-log(A/m/) (3.16) 



where A ~ 1/L is the UV cutoff of the effective GN model ( |3.4| ). Note that the dynamically 
generated mass scale 

m f ~ A e ~ 27T/X9n (3.17) 
is much smaller than the cutoff A for small GN coupling. 
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To complete the discussion, we need to show that the resulting f(x) is slowly varying 
on the scale L (to justify ( |3.13| )). Fourier transforming ( |3.14| ) we find the following behavior. 
For i-Ci, 

f{x) ~ /(0) = m f . (3.18) 

For L C z < 

while for a; > 1/m/, / goes exponentially to zero. 

Although f(x) for x > L is not constant, it is very slowly varying at weak coupling. 
Indeed, 



/(0) - f(x)\ x ^ L _ ln(xm f ) _ Xgn 



x 



/(0) " * ln(Lm/) " 2tt ln L • (3 ' 20) 

The logarithmic variation (|3.19| ), (|3.20| ) becomes important only on scales for which 
In £ ~ 1, i.e. x ~ 1/m/. We conclude that /(x) ( |3.11| ) varies on the scale 1/m/ ^> L. 
The same conclusion can be reached by examining the momentum space expression 
( p.!4j) . The Fourier transform f(k) contains two momentum scales, m/ and 1/L. However, 
most of its support is on the scale m/. At k = 1/L it has been reduced by a factor e _2L / A 
from its value at the origin. This means that the scale of variation of f(x) is 1/m/. 
The solution for / gives the following position dependent condensate: 

fat (x) . ,„(0)> = N c m f | |<a • (3.21) 

This is nothing but the massive propagator for a fermion of mass m/. Thus, we con- 
clude that the intersecting D-brane model generates dynamically a mass m/ ( |3.17|) for 
the fermions living at the intersections of the color and flavor branes. This mass is non- 
perturbative in the Gross-Neveu coupling ( jOl) . Comparing ( ^6|) to (|3.17|) we see that is 
exhibits the same dependence on the coupling as in the field theoretic analysis of 0. 

One can take a decoupling limit in which the physics of the brane configuration reduces 
precisely to the Gross-Neveu model. To do that one takes X/L — > and focuses on energy 
scales of order m / . This sends the ratio of the mass of the fermions m / ( [3.171) and the U V 
cutoff scale 1/L to zero, and leads to a model that is precisely equivalent to that of 0. 
The model with finite X/L can be thought of as a version of Gross-Neveu with a finite UV 
cutoff. 

In the above discussion we took the renormalization scale of the theory to be the 
UV cutoff A ~ 1/L. X gn ( p.5| ) is the coupling at that scale. One can choose some other 
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renormalization scale \i < A and define the coupling X gn at that scale. The renormalization 
group guarantees that Green functions expressed in terms of the physical fermion mass are 
insensitive to such changes. 

The Gross-Neveu model ( |3.4| ), which describes the dynamics of and q^ at weak 
coupling, was studied extensively in the past, and much is known about it. We next review 
some of its pertinent features. We first address the question of its extreme low-energy limit. 
This can be analyzed as follows (see e.g. The CFT of free massless fermions q^, q^ 

can be bosonized in terms of a WZW model for the group 

SU(N c ) Nf x SU(N f ) Nc x U(l) . (3.22) 

The subscripts in (|3.22| ) denote the levels of the corresponding current algebras. One 
can write the left-moving currents corresponding to ( |3.22| ) in terms of the fermions q^ as 
follows: 

SU(N c ) Nf : J1 = (ql)f(q L y p ~ j^SfJ ; (3 . 23) 
SU(N f ) Nc : Jj = (g£)$W a - j^Sp . 

Here a,f3 = 1,---,N C are color indices and i,j = l,---,Nf are flavor ones. A similar 
decomposition holds for the right- moving fermions. The currents on each line of ( |3.23| ) 



satisfy standard affine Lie algebra OPE's, and commute with those on the other two lines. 
The U(l) current on the first line describes a massless free field. The rest of the low-energy 
theory ( |3.22| ) is an interacting WZW model (for generic Nf, N c ). 

The four- Fermi interaction ( |3.4j) can be written in terms of the currents (|3.23|) , as a 



linear combination of an abelian Thirring interaction for the U(l) factor, and a non-abelian 
Thirring interaction corresponding to SU(N C ). The abelian Thirring interaction 

Ci = XiJJ (3.24) 

corresponds to changing the radius of the decoupled scalar field associated with the U(l) 
factor in (|3.22| ) . Thus, its effects are easy to analyze. The non-abelian Thirring interaction 
associated with SU(N C ), 

£ 2 = A 2 Jpi , (3.25) 

leads to non-trivial long distance dynamics. The coupling A 2 is asymptotically free. It 
grows in the infrared, where the SU(N c )N f degrees of freedom become massive and de- 
couple. The extreme infrared theory is a coset model. It is obtained from the UV theory, 
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which contains the current algebra ( |3.23j ), by gauging the SU(N c )]y f subalgebra. This 
eliminates the corresponding factor in ( |3.22| ) and leads to a WZW model for 



SU(N f ) Nc x U(l) . (3.26) 

One can think of this model as describing the infrared fluctuations of the order parameters 
associated with chiral symmetry breaking. 

As discussed in []6|,[7| these fluctuations effectively restore the chiral symmetry at finite 
N c . For example, the two point function 

(ql(x) ■ q R (x)J R (y) ■ q L (y)) (3.27) 

does not in fact approach a constant at large separation, as would be expected from ( |3.21|) , 
but rather decays like \x — Thus, at infinite iV c the symmetry is broken, but at 

finite N c it is restored at sufficiently long distances. As iV c increases, the range of distances 
for which the symmetry appears to be broken increases exponentially with iV c , and the 
large N c analysis becomes a better and better approximation to the dynamics. 

It will be useful below to determine the location of the massless modes ( |3.26| ) in the 
extra dimensions. This can be done following || by integrating out the fermions qL, qR 
and studying the normalizable modes of the resulting theory of gauge fields. To determine 
the fate of the SU (Nf) factori in the WZW model ( ^.26| ) it is enough to consider a single 
intersection, say that of the N c D4-branes and Nf D6-branes. Since we are interested 
in the low-energy dynamics, we can restrict to s-waves on the four-spheres in the (56789) 
directions and write an effective three-dimensional Lagrangian for the three-dimensional 
SU(Nf) gauge field with components! (A + , A_, A u ). That Lagrangian contains three 
types of terms, 

£ -^kin -£ferm -^inflow • 

(3.28) 



The first is the standard kinetic term, which has the form 



Adn = \ w 4 Tr 

97 



(3.29) 



The second comes from integrating out the fermions at the intersection. It is localized at 
u = and has the form 

£ ferm = 5(u)N c T(A(0)) . (3.30) 



2 The U(l) factor in (3.26) can also be treated as in 

3 u is the radial direction in (56789). 
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Here T(A) is the chiral Polyakov-Wiegmann action || corresponding to SU(Nf), and A(0) 
stands for the two-dimensional gauge field at the intersection, A±(x + , x~, 0). 



The third term in (|3.28|) is a Chern-Simons term, which arises as follows. On the 
worldvolume of the L>6-branes there is a coupling of the SU(Nf) gauge field A to the RR 
four- form field strength^ 

^2 J H 4 Au 3 (A) (3.31) 
where uj 3 {A) is the Chern-Simons form for SU(Nf), 

uj 3 (A) = Tr ^4 AF + ^A A A A A^j . (3.32) 

One way of deducing the presence of the term ( |3.31| ) and its coefficient is to require the 
cancellation of the SU(Nf) anomaly of the fermions living on the DA — D6 intersection by 
anomaly inflow flCfl . 

The N c D4-branes ( |2.1| ) are localized in the directions (56789) and are magnetically 
charged under H±. Therefore, there are N c units of H4 flux going through the four-sphere 
at fixed u, 

#4 = 2ixN c . (3.33) 



l S 4 

Integrating (|3.31[) over the sphere leads to a Chern-Simons term in the Lagrangian ( |3.28|) : 



^inflow = ^3 (A) . (3.34) 



The resulting Lagrangian ( 3.28 ) is very similar to one that was recently studied in ||. An 
analysis similar to that of Q shows that the SU(Nf) degrees of freedom, which classically 
live at the intersection u = 0, are displaced to u ~ . This displacement is small for weak 
coupling, but it has interesting implications at strong coupling, as in 0. 

Note that while in the above discussion we analyzed the position of the chiral SU(Nf) 
modes associated with a given intersection, the answer is the same for the system we are 
interested in, which has two intersections. As explained above, in that system there is a 
Thirring interaction for the SU(N C ) and U(l) factors in (|3.22| ). Since the SU(Nf) currents 
commute with these perturbations, they are not influenced by them. This decoupling 
between the massless sector and a non-trivial massive sector of a two-dimensional field 
theory is a rather general phenomenon; see e.g. [Tl\ for a discussion. 



We normalize the RR forms and Chern-Simons couplings as in 
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So far we focused on the properties of the massless modes associated with the brane 
configuration of section 2. The massive spectrum of the non-abelian Thirring model ( 3.251 ) 
is also known, since the model is integrable [p~2| - p~5[1 . The masses of single particle states 
are given by 



sin 



to/ 



sin 



(3.35) 



The integer n runs over the range n = 1, 2, • • • , N c — 1. The n = 1 states are the original 
fermions. States with higher n can be thought of as bound states of n fermions. 



4. Strong coupling analysis 

Bringing the L>6-branes closer or increasing g s N c increases X/L and takes us to the 
strong coupling regime in which the approximations used in the previous section break 
down. For L < A we can employ a dual description of the model by studying the DBI 
action for L>6-branes in the near-horizon geometry of the iV c L>4-branes. The discussion 
follows closely that of 

4-1. Chiral Symmetry Breaking 

The near-horizon metric and dilaton of the D4-branes are given by 

3 3 

ds 2 = (%)' [ n ^dx"dx- - (dx 1 )*] - (^) 2 (dU'< + f/ 2 dS]|) , 

where we set a' = 1, O4 labels the angular directions in (56789), and the parameter R is 
given by 

R s = 7vg s N c = ^-N c = nX . (4.2) 

Recall that the (4 + l)-dimensional 't Hooft coupling A has units of length. There is also 
a non-zero RR four- form field strength H4 with N c units of flux around the 5" 4 in Q4.1|) . 

The L>6-branei is stretched in the directions (01), wraps the four-sphere labeled by 
O4, and forms a curve U = U(x 4 ) in the (U : x 4 ) plane. The induced metric on it is 



3 . 3 r 



fu 



-3 



k/x 4 ) 2 -(^) 2 U 2 dQ\ 



3 



(4.3) 



We write all the formulae for Nf = 1, but it is easy to generalize them to larger Nf. 
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where U' = dU/dx 4 . The curve U(x 4 ) is determined by solving the equations of motion 
that follow from the DBI action for the sixbrane 



S D 6 = ~r Q V l+1 V,Ri J dx 4 U^ U" . 



(4.4) 



Vx+i is the volume of (1 + l)-dimensional Minkowski spacetime, V4 is the volume of a unit 
S 4 , and tq = T%jg s is the tension of a D6-brane for constant dilaton = g s . 

Since the Lagrangian does not depend explicitly on x 4 , there is a first integral given 

by 

If? 5 

U 2 (4.5) 



i + (S) u' 2 

where Uq is the value of U where U' = 0. We are looking for solutions which asymptotically 

L 
2 



approach a D6 and D6-brane a distance L apart, i.e. U(x 4 — > ±4) — > 00. Uq is the smallest 



value of U along the brane. 

There are two types of such solutions. The first is a straight brane and anti-brane 
stretched in U and located at x 4 = ±4, corresponding! to infinite U' and Uq = in (|4.5|) . 
This solution preserves the chiral SU(Nf)i, x SU(Nf)ji symmetry. 

The second is a symmetric U-shaped curve in the (U, x 4 ) plane. Integrating (|4.5|) , it 
is given by 



u du IRi 



^(U) = ±l -r = ± 



3 1 \ „ ( Vl 3 1 



5 ' 2 J V t/ 5 ' 5 ' 2 



(4.6) 



where -8(2; a, 6) is the incomplete Beta function, with S(0; a, b) = and 5(1; a, 6) = £?(a, 6) 
the usual Euler Beta function. The brane separation is 



2 Ri „ /3 1 



2x 4 (oo) = -— S -,- . (4.7) 



5rr^ V5 2 

This solution can be thought of as a D6-brane and a D6-brane connected by a wormhole, 
whose radius at its narrowest point is Uq. The fact that what looks asymptotically like D6 



6 Since the curvature diverges as U — > 0, the small U behavior of this solution is not reliably 
described by supergravity. This is not going to matter for our considerations that will mainly 
involve the large U region. 
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and D6-branes is in fact part of a single connected curved brane implies that the chiral 
U(Nf)L x U(Nf)n symmetry is broken down to the diagonal U(Nf). 

To find the ground state of the system one has to compare the energies of the straight 
and curved solutions for fixed L. The energy difference AE = -^straight ~ -^curved is pro- 
portional to 



AE ~ / du(u-Q)+ duu 

JUo 



u° 



\ulB (-\\) ^0.139 U%, 



(4.8) 

so the curved configuration is preferred and chiral symmetry is broken. We see that the 
system is in the same phase for weak and strong coupling. 

The regime of validity of the supergravity approximation is the same as in j2|,[[] , A ^> L 



or large effective £ t Hooft coupling, g s (Uo)N c ^> 1, [ 16|j 



4-2. Spectrum 

The wormhole connecting the DQ and D6-branes in the vacuum configuration de- 
scribed in the previous subsection provides a mass to the fermions qL,QR- These corre- 
spond at strong coupling to fundamental strings that stretch from the bottom of the curved 
flavor branes down to U = 0.@ Using the metric ( |4. 1| ) one finds that their mass is equal 
to Uo/2Tva' = Uq/2tt, which from ( \i.7\) is proportional to X/L 2 . This should be compared 
to the weak coupling regime, where we saw that the dynamically generated mass of the 
fermions ( |3.17| ) was given by rrif ~ (1/L)e -L / A . Defining 

nifL = A(X/L) , (4.9) 

A(x) behaves as 

A(x) ~ / e " 1/X ' x^O; (410) 
X, x — > oo. 



7 The fact that the lowest lying state of such a string is a fermion can be understood as follows. 
The background Q4.1D is close to flat except very close to U = 0. The open string in question ends 
on one side on a brane which is stretched in the directions (01), x A and an S A of radius Uo- The 
other end of the string is at U = 0, where the curvature is large, but we can replace this region by 
boundary conditions corresponding to an open string ending on D4-branes extended in (01234). 
Overall, the string has six Dirichlet-Neumann directions (the directions (23) and the S 4 ), so the 
lowest lying state is a (massive) fermion. 
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For x ~ 1 both our approximations break down, but it is natural to expect A(x) to be 
monotonic in x, so that the dynamically generated mass is a monotonic function of the 
coupling. 

In addition to the long strings discussed above the spectrum includes light mesons 
which correspond to normalizable fluctuations of the scalars, gauge field and fermions on 
the D6-brane. To study them we need to expand the DBI action around the background 
configuration (|4.6|) . In the D4 — D8 — D8 system it was shown in [|] that fluctuations of 
the gauge field give rise to massless Nambu-Goldstone (NG) bosons and massive vector 
mesons, while scalar fluctuations lead to a spectrum of massive scalar mesons. We will 
see that in our model gauge field fluctuations do not give rise to normalizable NG modes. 
Also, while the kinetic term of the gauge field gives rise to a massless vector, it acquires a 
mass by mixing with other modes via the Chern-Simons interaction (|3.31| ). 

We start with the action for the probe D6-brane including the worldvolume gauge 
field: 



~>D6 



-T 6 J d 7 xe-*^-det(g AB + 2nF AB ) + ^ J H 4 Au 3 (A) . (4.11) 



Here is the RR field sourced by the color L>4-branes (as in (|3.31|) ) and 1V3 is the 



Chern-Simons three- form ( |3.32| ). In the four-dimensional model of |§ there is an analo- 
gous coupling to the Chern-Simons five-form, which plays an important role in anomaly 
cancellation, but since it is cubic in the gauge field, it does not affect the mass spectrum. 
In our case, the Chern-Simons term contains terms quadratic in the gauge field and must 
be included in the analysis. 

We focus on SO(5) singlet gauge field fluctuations that are independent of the S 4 
coordinates. Denoting the remaining components of the gauge field by A m , m = 0, 1,U 
and expanding the action to quadratic order we find 

5 (quad) = _ T6 y 4(27r )2 J d^dUe-* ^fetii- A F mn F mn + ^ J d 2 xdUe mnp A m F np . 

(4.12) 

In this expression detg is the determinant of the induced metric on the D6-brane (|4.3|) , and 
the indices m, n are raised and lowered with the m, n components of this metric. Using 
the solution ( |4.6| ) we can rewrite the action ( [4.12| ) as 



^(quad) = _ T6 y 4(27r) 2 J d * xdU \\ G (JJ) F ^ F ^ + Ih{U)F^uF^u 

J d 2 xdUe^ (AuF^ + 2A„F uU ) . 



47T 



(4.13) 
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In ( 4.13 ) the indices //, v = 0, 1 are raised and lowered with the flat metric r\p V and all U 
dependence is contained in the functions 



G(U) 



9sU 2 



H{U) = ~~~~tt , (4.14) 



7(17) 



U 



5 



We now separate variables and expand the gauge field in modes in the U direction as 
follows (with normalization conditions to be determined shortly) 



A^x\U) = Y,B { ;\x v )MU), 

n 



(4.15) 



Defining B$ = d^B^ — d v B^ this leads to the action 



>(m)/ii/ 



S = -T 6 V 4 (2n) 2 I d 2 xdU^G{U)^J- A B^B 

n,m 

+ ^H{U) (<V (n Vn - Bjfiduipn) - B^du^ 

n,m 



n,m 



(4.16) 

We see from the first line of ( |4.16|) that the kinetic term for is canonical provided 
that the wavefunctions i^ m {U) satisfy the normalization condition 

T 6 V 4 (2n) 2 J dUG{U)^ n {U)^ m {U) = 5 n , m . (4.17) 

Terms on the second line of (|4.16|) give diagonal and canonically normalized mass terms 
for Bp with mass squared provided that we also have 

T 6 V 4 (2n) 2 J dUH(U)du^ n (U)dutj} m (U) = m 2 J n , m . (4.18) 



Equations (|4.17| ) and ( |4.18| ) are compatible if ip n satisfies the eigenvalue equation 

-^du [H(U)duMU)} = -mlMU) • (4.19) 
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Finally, in order to have canonically normalized kinetic energy terms for the we require 

T 6 V 4 (2n) 2 J dUH(U)<p n (U)<Pm(U) = S n , m . (4.20) 

This can be satisfied by choosing <fi n = duil>n/\iTi>n\ for modes with m n 7^ 0. For m n = 0, 
a formal extrapolation of the above analysis leads to a zero mode, <fio ~ d u ipQ ~ 1/H(U). 
In the D4 — D8 — D8 system, an analogous mode was shown in to be the massless 
Nambu-Goldstone boson associated with the breaking of U (Nf) l xU (Nf)n to U (iV/)di a g- 
In our case, in contrast to 0, this putative zero mode is neither orthogonal to the other 
<fi n , nor normalizable. The latter follows from the fact that H(U) ~ U at large U so that 
the integral 

' dUH(U)MU)MU) ~ / TF777V ( 4 -21) 



H(U) 



diverges logarithmically at large U. Thus, there is no normalizable zero mode and hence no 
massless NG mode in the spectrum. This seems odd, since we have established above that 
the symmetry is broken and one expects this to lead to a NG boson on general grounds. 
We will discuss the resolution of this puzzle below. 

Another puzzling feature of the spectrum is the existence of a normalizable zero mode 
for the gauge field, corresponding to ^0 = const in (|4.15|) . It satisfies (|4.19| ) with m n = 
and from Q4.17 ) is normalizable since G(U) ~ 1/U 2 at large U. This seems to suggest that 
the U(Nf) symmetry on the D6-branes is gauged in 1 + 1 dimensions, which is unexpected 
from the weak coupling point of view. We will see below that this mode is not an eigenstate 
of the full quadratic action (^4.16|) . 

To summarize, in the expansion (|4.15| ) there is a normalizable zero mode i/jq in the 
expansion of A^, but the expansion of A\j involves a sum only over non-zero modes 
for which <fi n = du^ n /\ m n\- We thus see that Ajj is pure gauge, Ajj = djjK with 
A = J2 n j^~\ 7i; ^ n \ x )' l l } niU). Making a gauge transformation by A to set Ajj = and 
substituting Q4.15 ) into the action ( 4.13 ) leaves us with the two-dimensional action 



Sdb = - d x 



1 B (0) B ^(0) + y l _ B (n) B ^{n) + 1 2 B (n) B 
n>l V 

d 2 xJ2^M 2 mn B^B^ 



(4.22) 



where 



M 2 nm = -M 2 mn = ^ j dUrjUUh r,jU). 



(4.23) 
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It is clear from Q4.22 ) that the fields B^ are not mass eigenstates because of the off- 



diagonal mixing. It is difficult to compute the mass eigenstates explicitly since to do so we 
would have to diagonalize the infinite-dimensional matrix with diagonal components m 2 
and off diagonal components .M 2 m . However, there is no reason why the determinant of 
this mass matrix should vanish, so this process presumably leaves us with mass eigenstates 
all with non-zero masses. In particular, the massless gauge field mixes with the 

massive modes via (f4.23| ) with n = 

AT f 

M 2 0m oc / dUduip M (U) . (4.24) 
Zn J 

It is easy to check that normalizable solutions of ( |4 . 1 9| ) approach a non-zero constant as 
U — > oo, so the integral in (4.24), and hence the mixing of B^ with the other modes, is 



non-zero. 

One general observation that can be made about the spectrum is the following. It was 
recently emphasized in [T7J that the holographic distance-energy relation of |H| coupled 



with the consistency of holography for open and closed string modes implies a universal 
scaling of meson masses in the strongly coupled, supergravity regime: 

M meson (4.25) 

where m/ is the fermion mass and (7 2 ff is the effective 't Hooft coupling evaluated at the 
scale m/ . In our model it is given by 

9es( m f) - gt N cmf ~ Am/ . (4.26) 



Using the fact that m/ ~ A/L 2 , (|4.25|) leads to the prediction that all meson masses are 



proportional to M meS on ~ 1/L and are thus deeply bound since M meS on 'Cm/. To see 
that this is in fact true, note that the mass parameters m 2 in ( |4.22| ) are determined by 
the eigenvalue equation (|4.19| ). Writing it in terms of U = U/Uq gives 



U 7 " 5 % C/7" 3 %^n = -Trm^n . (4.27) 



Since all the dependence on the parameters A and L is now contained in the factor R 3 /Uo ~ 
L 2 , we see that the eigenvalues m 2 are all proportional to 1/L 2 and independent of the 't 
Hooft coupling A. 
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The mixing parameters -M^ n ([4.23 ) also have this property. The normalization con- 
dition for ip n ( f4.17| ) can be written as 

%i^3 m = Km (4.28) 

u 



where 



*• = T7m ,K (4 ' 29) 



" c 

for c a constant, independent of L and A. Therefore, n ( |4.23| ) also scale like 1/L 2 . 

To summarize, fluctuations of the gauge field lead to a spectrum of massive vector 
fields in 1 + 1 dimensions. The masses all scale like 1/L and are determined by solving the 
eigenvalue equation Q4.19|) and then diagonalizing the quadratic form in Q4.22| ). A similar 
analysis leads to a spectrum of massive scalar and fermion fields with masses that scale 
like 1/L. 

In contrast to the model of Sakai and Sugimoto 0, we do not find any massless NG 
bosons in the spectrum, in spite of the apparent breaking of U(Nf)i, x U(Nf)n to the 
diagonal U(Nf). The resolution of this puzzle was pointed out in a closely related context 
in ||. We saw in section 3 (after eq. (|3.34|) ) that in the weak coupling limit, the NG 
bosons are displaced from the intersection of the branes by an amount of order A 1 / 3 in 
the U direction. Extrapolating this to strong coupling we find that the U(Nf) degrees of 
freedom live in the transition region between the near-horizon geometry of the L>4-branes 
and the asymptotically flat space far from the branes. Therefore, they are not visible in 
an analysis that restricts attention to the near-horizon geometry. Note that this kind of 
behavior is only possible in two dimensions where the NG bosons are decoupled from the 
massive physics. In higher dimensions similar behavior is familiar from the analysis of the 
center of mass modes in brane systems, which decouple in the near-horizon limit. 

The scalar, fermion and vector mesons obtained by studying small fluctuations of the 
D-branes are quite different from the bound states of the GN model whose masses are 
given by (|3.35|) . In particular, as is clear from our analysis, at weak coupling the color 
degrees of freedom living on the L>4-branes are essentially non-dynamical, and their main 
effect is to produce the small attractive force between the fermions q^, qR- AH the states 
( P-35|) can then be thought of as bound states of these fermions. 

On the other hand, at strong coupling, the degrees of freedom living on the L>4-branes 
(or, more properly, compactified M5-branes) are strongly interacting, and the states with 
masses of order 1/L are best thought of as containing the fermions qL,R as well as the 
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color degrees of freedom. They have the qualitative form q^<& n q, where q stands for the 
fermions living at the intersection, and $ stands for all the bosonic and fermionic degrees 
of freedom on the D4-branes. Therefore, these states do not have a direct analog at weak 
coupling. It is interesting that they are much lighter than the fermions themselves, like 
Nambu-Goldstone bosons of a spontaneously broken symmetry. 

To study the analogs of the GN states (|3.35| ) at strong coupling, one needs to consider 
long fundamental strings stretched between the curved L>6-branes and U = 0. As men- 
tioned above, a single string of this sort is a massive fermion. A few such strings might 
form bound states, like their weak coupling counterparts. This is an interesting problem 
for future study. 



5. Finite temperature 

At zero temperature, the intersecting brane model of sections 2-4 exhibits dynamical 
chiral symmetry breaking both at weak coupling (discussed in section 3), where it reduces 
to the Gross-Neveu model, and at strong coupling, where it is well described by probe D6- 
brane dynamics in the near-horizon geometry of L>4-branes (section 4). In this section we 
will study this model at finite temperature, and in particular analyze the finite temperature 
phase transition at which chiral symmetry is restored. In the following two subsections we 
do that in turn for weak and strong coupling. In the last subsection we briefly address the 
interpolation between the two regimes. 

5.1. The Gross-Neveu model at finite temperature 

The finite temperature properties of the Gross-Neveu model were studied in [|T9| - pT | 



and many subsequent papers. Here we briefly review the results, which are relevant for our 
intersecting D-brane system in the limit \ gn — ► 0, E ~ mj discussed in section 3. In this 
limit the D-brane system is described by the local Lagrangian ( |3.4| ), which can be written 
in a way similar to 



S e & = j d 2 x (^iq j L a^d fl q L + iq R a fl d fl q R + aq R - q L + aq\- q R - J^ (J(J ^ ■ 

The expectation value of a is the mass of the fermions qL,<lR- To calculate it, it is con- 
venient to integrate out the fermions and minimize the resulting effective potential for 
a. 
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For zero temperature, the effective potential takes the form 

1 



eff 



A 



-aa 



"a 



d k , / aa 



(5.2) 



Stationary points of V e Q satisfy 

SVeS 

5a 



a 



A 



d 2 k 



a 



{2n) 2 k 2 + aa 



(5.3) 



which is equivalent to the gap equation ( |3.15|) discussed in section 3 in the local, GN, 
limit, a = is a local maximum of the potential Q5.2p; the minimum corresponds to the 



non-trivial solution of (|5.3| ) where chiral symmetry is broken. 

To study chiral symmetry breaking at finite temperature, we need to evaluate the 
expectation value 

(*> = (5-4) 



where H is the Hamiltonian and 1//3 = T the temperature. This can be done by studying 
the theory corresponding to ( |5.1|) on IR x S , with the Euclidean time living on a circle 
of circumference f3 with anti-periodic boundary conditions for the fermions ql, qn. This 
amounts to replacing in ( |5.2| ), (Of), 



dko 



- r 



(5.5) 



with ko = 2nr/P the momentum in the Euclidean time direction. Making this replacement 
in ( |5.2j ) and minimizing the effective potential with respect to a gives the leading contribu- 
tion to the free energy F(f3). The value of a at the minimum gives the expectation value 



The stationarity condition of V e s is now 



SVeff 

5a 



a 



A 



gn 



1 _ ~^3 n 
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E 



dp 



A 



2tt 



2nr 

S3 



+ p 2 + \a\ 



(5.6) 



Following one can replace the sum over r by a contour integral over energy, 

de 1 



- E 



2ni eP e + 1 



(5.7) 



22 



where the contour surrounds the poles at f3e = 2-nri. Deforming the contour and picking 
up the residues of the poles of the r.h.s. of (|5.6|) at 



±e p = ±VP + W 



(5J 



gives 



SVeS 

5a 



a 



A 



4tt 7_ a e 



dp 

p 



— tanh ( -(3e p 



(5.9) 



The stationary points are a = and the solution of 



2tt 



dp 



tanh I —j3e 



As (3 — > oo we can approximate tanh ^/9e p ~ 1, and ( |5.10|) (with cr(/3 — > oo) 



1 = v r^p^ v ln ,2A 



2tt 



27T 771/ 



(5.10) 
m/) reads 
(5.11) 



Thus, the mass of the fermions at zero temperature is given by 

m/ = 2Ae A f n , 



(5.12) 



which agrees with (|3.17| ) after allowing for a different definition of the cutoff A. 

For low temperature, i.e. large but finite (3, a = is still a local maximum. Indeed, 
for a = 0, e p = \p\ and the integral in ( |5.9|) is given by pl| : 



f A dv 1 

/ — tanh(-/3p) ~ ln(1.14/3A) 
Jo P 2 



(5.13) 



where we assumed that f3A 3> 1, as appropriate in the local GN limit. Thus, one has 

?2 



5 2 V eS (a = 0) _ 1 ^ A gn 



A 



3" 



ln(1.14/3A) ] = -^-ln(0.57/3m f ) , (5.14) 

27T y 27T 



where in the last step we used ( |5.11| ). We see that a = is a maximum of V e s for (3 > (3 C 
and a minimum for (3 < (3 C . The critical temperature is 



1 



0.57m/ 



(5.15) 



Below that temperature, the minimum of the effective potential is at a non-zero value of 
a, a = rrif which corresponds to a non-trivial solution of ( |5.10| ) with chiral symmetry 
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broken. As T — > T c , mf((3) — > 0; for T > T c the only solution of ( |5.10| ) is a = 0, and 
chiral symmetry is restored. The phase transition at T = T c is second order since the order 
parameter m/(/3) is continuous at the transition. 

It should be noted that, as in section 3, the discussion above is strictly speaking valid 
only for infinite N c . For finite N c the order parameter ( |5.4|) is not really constant in space 
as was assumed in our analysis. Averaging over sufficiently large distances, this expectation 
value vanishes for all T. However the domains over which it is roughly constant grow with 
N c . Thus, the large N c analysis is a useful guide to the behavior of the system for finite 



N c as well (see [|I9| - ^I| for further discussion). 



5.2. Strong coupling analysis at finite temperature 

At strong coupling the dynamics of the intersecting branes is described by studying 
probe D6-branes in the near-horizon geometry of the N c D4-branes. To study the system 
at finite temperature we compactify Euclidean time on a circle of circumference (3, with 
anti-periodic boundary conditions for the fermions. The resulting metric and dilaton are 
given by 



3 3 



ds 2 = QQ 5 (f(u)(d X r + (dxr + (dxr) + (|) " 5 ^ 2 + ^ 



e = 9s I ■= 



with 



(5.16) 

f(U) = 1 - ^| • (5.17) 

This geometry is a continuation to Euclidean space of the non-extremal D4-brane back- 
ground. The non-extremality parameter Ut is related to the temperature via 

(3= -. (5.18) 

3U% 



For zero temperature one has Ut = 0, and the solution ( ^.16|) , ( |5.17|) reduces to (|4.1|) 



To calculate the free energy of the intersecting brane theory we need to solve for the 
shape of the D6-branes in the geometry ( |5.16| ). This can be done similarly to the closely 
related DA — D8 — D8 case which was studied in |22|,[23 . 
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The DBI action for the sixbranes is 



S 



D6 



-TeVi 



+1 



I 



dx 4 U^lf(U)+(^) If' 2 



(5.19) 



The generalization of ( |4.5[ ) to non-zero temperature is 



U' 2 
f(U) 



(5.20) 



Uq is the lowest value of U reached by the brane. 

We need to solve ( |5.2U| ) for the shape of the branes U(x 4 ), as in section 4. There are 
again two possible solutions. In one the D6 and D6-branes wrap the smooth space labeled 
by (U, x ), at two different values of x A a distance L apart. In this configuration the branes 
extend all the way to the Euclidean horizon; thus Uq = Ut- The chiral U(Nj)l X U(Nf)n 



symmetry ( |2.3|) is preserved. 

In the other solution, the D6 and L>6-branes are connected, with U = £7(x 4 ), and 
Uq > Ut- This solution breaks the chiral U(Nf)i, x U(Nf)R symmetry to the diagonal 
U(N f ). 

To see which of the solutions has lower free energy we need to compare the values of 
the corresponding effective potentials 



AE = E. 



straight 



E, 



curved 



U 



dU(U-0)+ / dUU 



U T 
1 

VT 



U 



1 



1 



U$f(Uo) X 
U 5 f(U) J 



dy(y-0)+ / dyy 



1 - 



/(I) 

y 5 f(y) 



' 2 



(5.21) 

where in the last line we introduced the dimensionless coordinate y = U/Uq. 

The numerical result for AE(T) is plotted in figure 2. For low temperature, the curved 
solution is energetically preferred. AE vanishes at a critical temperature corresponding to 
y T = U t /Uq ~ 0.613, above which the symmetric solution becomes energetically preferred. 
Thus, the system exhibits a first order phase transition at a temperature T c ~ 1/L. The 
constant of proportionality can be fixed as follows. In the low temperature phase, L, Ut 
and Uq are related by 



2^r / dy 



y 



(5.22) 
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At the critical temperature yr = y^p one finds L m 1.097 ^/2 , which implies 



3/2 



T c « 0.205/L . 



(5.23) 



The phase transition at T c is first order since the fermion mass jumps from a non-zero 
value to zero as T crosses T c . The mass right below the transition is equal to the energy 
of a string stretched from the curved D6-branes to the Euclidean horizon at U = Ut, 



m f 
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Fig. 2: Difference between the free energy of the straight and curved solutions as 
a function of temperature. 



5.3. Interpolating between weak and strong coupling 

In the last two subsections we discussed the finite temperature phase transition in the 
limits X/L — ► and X/L — > oo. It would be interesting to understand the interpolation 
between the two regimes. 

The transition temperature, T c , is given in the two limits by ( |5.15| ) and ( |5.23| ), respec- 
tively. A natural interpolation to all values of the coupling is 



^ = H(m f L) . 
m f 



(5.25) 



Here rrifL is a measure of the coupling X/L. It goes to zero in the GN limit and to infinity 
in the supergravity one (see (|4.9|), (|4.10|) ). H is a function to be determined. Eq. (|5.15Q 
implies that 

fT(0) = 0.57 (5.26) 
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while ( |5.23|) implies that 

205 

H(x -> oo) ~ — - . (5.27) 

x 

It is natural to expect that H(x) varies smoothly and monotonically with x. 

Another interesting observable is the order parameter right below the transition, 

^§2 = B(m f L) . (5.28) 
m/(0) 

From the GN analysis we know that 

B(0) = , (5.29) 

i.e. the symmetry restoration transition is second order. We expect the transition to 
remain second order for sufficiently small X/L, assuming that the physics at finite X/L is 
continuously connected to that of the GN model.! The reason is that the effective potential 
( |5.9|) behaves near the origin like 

V eS = ai \a\ 2 + a 2 \a\ 4 + ■ ■ ■ (5.30) 



where a\ is given by ( |5.14[ ) and changes sign at T = T c and a 2 is strictly positive. Small 



X/L corrections to a± and a 2 cannot change the fact that the transition happens at the 
point where a± changes sign and is second order. 
At strong coupling we found 

B(oo)= 0.387. (5.31) 
We expect B to vanish for X/L below some critical value, and then smoothly interpolate 



between zero and (5.31) as X/L — > oo. 



6. Discussion 

In this paper we analyzed a model of intersecting D4 and D6-branes and antibranes 
that has the interesting property that the chiral symmetry of the original brane configu- 
ration ( |2.3[ ) is broken to the diagonal subgroup by the dynamics. The model contains an 
adjustable parameter, X/L, which governs the strength of the coupling between left and 

8 We thank O. Aharony for a discussion of this issue. 
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right-moving fermions living at different intersections. For weak coupling it reduces to the 
Gross-Neveu model ||, a vector model which is exactly solvable^ at large iV c . For strong 
coupling it can be analyzed by studying the DBI action for D6-branes in the near-horizon 
geometry of N c D4-branes and again breaks chiral symmetry. 

As one would expect, the strength of chiral symmetry breaking increases with the 
coupling. For X/L <C 1, the dynamically generated fermion mass mf ( |3.17| ) is much 
smaller than the natural scales in the problem, 1/L and 1/A. Indeed, m/L, m/A go to 
zero as the coupling goes to zero. For X/L ^> 1 the fermion mass is of order X/L 2 and 
is much larger than the other scales: m/L, m/A go to infinity as the coupling goes to 
infinity (see (|4.9|) , ( f4.10|) ). The strong interactions between the fermions and the degrees 
of freedom living on the D4-branes give rise to strongly bound mesons whose masses are 
much smaller than m/ (their binding energy is thus almost 100%). These bound states 
can be studied by expanding the DBI action of the sixbranes around the vacuum solution. 

Another manifestation of the dependence of the strength of chiral symmetry breaking 
on the coupling is the thermodynamic properties described in section 5. For strong cou- 
pling, the transition is strongly first order - the fermion mass right below the transition 
is of order m/ (see fl5.24D ), which as mentioned above is large in natural units. As the 
strength of the interaction decreases the transition becomes weaker and it is second order 
at weak coupling (as described in subsection 5.1). 

There are many possible directions for further investigation of this and related models. 
The Gross-Neveu model has a spectrum of meson excitations with masses given by ( |3.35| ) 
whose S-matrix is known due to integrability. In our setting, the GN model is obtained 
in the limit X/L — > described in section 3. It would be interesting to see whether the 
integrability persists for finite X/L, and if so what one can learn about the model by using 
it. In particular it would be interesting to determine the mass spectrum, S-matrix and the 
thermodynamic functions H ( |5.25| ) and B ( |5.28[ ) for all values of the coupling. 

Another issue involves the stability of the system beyond the large iV c limit. In the full 
IIA string theory the brane configuration of figure 1 is unstable. The D6 and L>6-branes 
attract each other by massless closed string exchange and will start moving towards each 
other if placed in this configuration. These gravitational effects turn off in the limit g s — > 0, 
and therefore can be neglected in the 't Hooft limit that we took. It is natural to ask what 
happens for finite iV c . We suspect that the system is still stable although this deserves 



In fact this model is integrable at finite iV c as well. 
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more careful study. On the supergravity side, taking the near-horizon limit means that the 
D6 and L>6-branes become infinitely separated near the boundary at U = oo. Therefore we 
expect the force per unit worldvolume to go to zero, while the mass per unit worldvolume 
remains constant. On the field theory side the dynamics is governed by the (2, 0) theory 
compactified on a circle. The D6-branes are interpreted as defects in this theory and we 
see no clear reason why the introduction of defects should result in a drastic instability of 
the theory. 

There are many natural generalizations of our model. The brane configuration of 
figure 1 contains a single stack of L>6-branes all at the same point in the IR 3 labeled by 
(234), and Nf antibranes at a distance L in IR 3 . As a consequence, the dynamics depends 
on a single parameter, X/L. One could generalize this construction by putting the D6- 
branes and antibranes at arbitrary points in IR 3 , and studying the dynamics as a function 
of their relative positions. At weak coupling this gives a generalization of the GN model, 
while at strong coupling it gives a generalization of the DBI dynamics of section 4. In 
both limits there is a rich phase structure with different patterns of symmetry breaking 
depending on the positions of the branes and antibranes in IR 3 . It would be interesting to 
analyze it. 

Another generalization involves making the SU(N C ) gauge fields dynamical. This can 
be done by compactifying the directions (234) on a torus of volume V234. The resulting 
model is very similar to the 't Hooft model of two-dimensional QCD ||24j| . The two- 
dimensional 't Hooft coupling, which has units of energy squared and sets the scale of the 
masses of mesons, is given by 

* = £ • «"> 

The dynamics depends on the hierarchy of scales and the strength of the various couplings. 
At weak coupling and large volume V234, the scale of confinement y/X~2 is much smaller 
than that of chiral symmetry breaking mj ( |3.17|) . Thus, at the scale of chiral symmetry 
breaking one can neglect the dynamics of the two-dimensional gauge field, and the mass 
generation of the fermions occurs as in section 3. At the much lower energy scale y/X~2 the 
effect of confinement kicks in, and the analysis of [23] is applicable (with fermion mass rrif). 
As the volume V234 decreases the physics changes, and for A2 3> mj the gauge dynamics 
becomes strong at energy scales at which the fermions are still massless. 

Thus, in this model confinement and chiral symmetry breaking are distinct phenomena 
and can occur at widely separated energies. This property was noted in a four-dimensional 
intersecting brane system in p| and is a general feature of such systems. 
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Other brane configurations with qualitatively similar dynamics are obtained by chang- 
ing the dimensionalities of the color and flavor branes and that of the intersection. This 
leads to many interesting models, including that of one that was considered recently 
in [25], and many others. We will discuss some features of such models in a companion 
paper |26| . 
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